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Abstract In this paper we are discussing various aspects of wavelet filters. While 
there are earlier studies of these filters as matrix valued functions in wavelets, in 
signal processing, and in systems, we here expand the framework. Motivated by 
applications, and by bringing to bear tools from reproducing kernel theory, we point 
out the role of non-positive definite Hermitian inner products (negative squares), 
for example Krein spaces, in the study of stability questions. We focus on the non- 
rational case, and establish new connections with the theory of generalized Schur 
functions and their associated reproducing kernel Pontryagin spaces, and the Cuntz 
relations. 
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1 Introduction 



Roughly speaking, systems whose inputs and outputs may be viewed as signals are 
called filters. Mathematically, filters are often presented as operator valued functions 
of a complex variable. In applications, filters are used in areas as (/) prediction, (ii) 
signal processing, (Hi) systems theory and (z'v) Lax-Phillips scattering theory ll55l . 
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There, one is faced with spectral theoretic questions which can be formulated and 
answered with the use of a suitable choice of an operator valued function defined 
on a domain in complex plane; in the case of scattering theory, the scattering oper- 
ator and the scattering matrix; in the other areas, the names used include polyphase 
matrix, see e.g., l43l |49l . We also mention that more recently, filters are used in 
(iv) multiresolution analysis in wavelets. We follow standard conventions regarding 
time-frequency duality, i.e., the correspondence between discrete time on one side 
and a complex frequency variable on the other. In the simplest cases, one passes 
from a time series to a generating function of a complex variable. These frequency 
response functions fall in various specific classes of functions of a complex variable; 
the particular function spaces in turn are dictated by applications. Again, motivated 
by applications, in our present study, we adopt a wider context for both sides of the 
duality divide. On the frequency side, we work with operator valued functions. This 
framework is relevant to a host of applications, and we believe of independent in- 
terest in operator theory. From the literature, we mention l57l . Il22l (see also ET\ ). 
and the papers referenced below. 

We here consider the set of C NxN - valued functions meromorphic in the open unit 
disk O0 and define two subsets of it: We shall denote by the family satisfying 
the symmetry 

W(e N z)=W(z)P N , (1) 
where En = e~®~ and denotes the permutation matrix 

p N= f Q lx(N-l) 1 

We shall also denote by °i/ lN the set of C A ' xA '-valued functions which take unitary 
values on the unit circle T. 

Classically wavelet filters, denoted by Wn, are characterized by rational functions 
satisfying both symmetries, i.e. 

W N = <% lN n<tf N . (3) 

In a previous paper, see j9j, we have provided an easy-to-compute characterization 
of Wn as both a set of rational functions, and in terms of state space realization. 

The aim of this work is to explore the possibility of extending the notion of wavelet 
filters, described in (|3). The functions considered still satisfy the symmetry in (HJ, 
but: 

• The functions are not necessarily rational or finite dimensional. 

• The functions are not necessarily unitary on the unit circle T. 



1 Classically, in the engineering literature, the functions are analytic, or more generally meromor- 
phic, outside the closed unit disk. The map z >-> 1/z relates the two settings. 

2 For rational functions, the term para-unitary is also used in the engineering literature. 
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• The functions are meromorphic (rather than analytic) in D. 

To explain our strategy, first recall the following: If W is a <C NxN - valued function 
which is rational and takes unitary values on the unit circle, the kernel 

i N -w(z)w( w y 

Kw(z,w)- — 

is positive definite in the open unit disk ID) if W has no poles there, or more gener- 
ally has a finite number of negative squares in D. See Definition 13 .41 below for the 
latter. In our approach, unitarity on the unit circle is replaced by the requirement 
that W is a generalized Schur function, in the sense that W is meromorphic in D and 
the associated kernel Kyy(z,w) has a finite number of negative squares there. This 
family includes in particular the case of matrix-valued rational functions which take 
contractive values on the unit circle. We will also consider the case where the values 
on the unit circle are, when defined, contractive with respect to indefinite metrics. 
These kernels are of the form 

j 2 -w( z )JiW( w y 

1 — zw* 

when W is C PlXpi -valued and analytic in a neighborhood of the origin, and where 
J\ and J2 are signature matrices, respectively in C Pl X,H and C P2Xpi , which have the 
same number of strictly negative eigenvalues: 

v_(7i) = v_(7 2 ), (5) 

and such that the kernel Kw has a finite number of negative squares. In |9) we stud- 
ied the realization of wavelet filters in the C ArxM -valued (with M >N) rational case. 
The above approach allows us to extend these results to the case where the filter is 
not necessarily rational and M may be smaller than N. Furthermore, the conditions 
in (2l of the function being analytic in the open unit disk, and taking coisometric 
values on the unit circle, are both relaxed (in particular, in the previous case, in ©, 
we had J\ = Im and J 2 = In)- 

The paper is organized as follows. Since we address different audiences, Sections 
2,3 and 4 are of a review nature. In Section [2] we give background on the use of 
filters in mathematics. We note that the more traditional framework in the literature 
has so far been unnecessarily restricted by two kinds of technical assumptions: (i) 
restricting to rational operator valued functions, and (/;) restricting the range of the 
operator valued functions considered. In Section[3]we address indefinite inner prod- 
uct spaces, and survey the theory of Pontryagin and Krein spaces. This overview 
allows us in Section [4] to describe a setting that expands both the above mentioned 
restrictions in (i) and (ii), namely the theory of generalized Schur functions. Our 
results in Sections 151 and l6l (Theorems 15 . 31 15 .41 and 16.41 ) deal with representations. 
We use these results in obtaining classifications, and decomposition theorems. In 
Section|71 we employ these theorems in the framework of wavelets. 
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2.1 Cuntz relations 

The Cuntz relations were realized by J. Cuntz in ll24) as generators of a simple 
purely infinite C* -algebra. Since then, they found many applications, and the related 
literature about Cuntz relations has flourished. Since Cuntz's paper (24), the study 
of their representations has mushroomed, and now makes up a big literature, see for 
example |[T9] [18j |20] [25] [J3] [37] [39l , and some of their applications |38ll40ll4T||42l 
[32), for example to fractals 1131") . 

In the initial framework, one is given a finite set Si , . . . ,Sjf of isometries with orthog- 
onal ranges adding up to the whole Hilbert space. Their representations play a role 
in a variety of applications, for example wavelets, and more generally multi-scale 
phenomena. The study of what are called non-type / C* -algebras was initiated in the 
pioneering work of Glimm l34l 135) and Dixmier 128) . This in turn was motivated 
by use of direct integrals in representation theory, both in the context of groups and 
C* -algebras. Direct integrals of representations are done practically with the use of 
Borel cross sections. Glimm proved that there are purely infinite C* -algebras which 
do not admit Borel cross sections as a parameter space for the set of equivalence 
classes of irreducible representations; the Cuntz algebra(s) On is the best known 
examples, 124) . Nonetheless, it was proved in |fT8) that there are families of equiv- 
alence classes of representations of On indexed by wavelet filters, the latter in turn 
being indexed by infinite-dimensional groups. 

One illustration of the need for expanding the framework of On from Hilbert space 
to the case of Krein spaces is illustrated by applications to scattering theory for the 
automorphic wave equation 154) . The initial study was restricted to the case when 
the operators 5, act on Hilbert space, and when they act isometrically. However, 
since then, there has been a need for generalizing the Cuntz relations. It was noted 
in |[T9) that the isometric case adapts well to the restricted framework of orthogonal 
wavelet families 1126) . Nonetheless, applications to engineering dictate much wider 
families, such as wavelet frames. 

In this work we extend what is known in the literature in a number of different di- 
rections, including to the case of Pontryagin spaces. We obtain Cuntz relations for 
isometries between certain reproducing kernel Pontryagin spaces of analytic func- 
tions. 



2.2 Wavelet filters 

In electrical engineering terminology, systems whose inputs and outputs may be 
viewed as signals are called filters. By filter, we here mean functions W (z) defined 
on the disk in the complex plane and taking operator values, i.e., linear operators 
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mapping between suitable spaces. 

While filters (in the sense of systems and signal processing) have already been used 
with success in analysis of wavelets, so far some powerful tools from systems the- 
ory have not yet been brought to bear on wavelet filters. The traditional restriction 
placed on these functions W(z) is that they are rational, and take values in the uni- 
tary group when z is restricted to have modulus 1 . In models from systems theory, 
the complex variable z plays the role of complex frequency. A reason for the recent 
success of wavelet algorithms is a coming together of tools from engineering and 
harmonic analysis. While wavelets now enter into a multitude of applications from 
analysis and probability, it was the incorporation of ideas from signal processing 
that offered new and easy-to-use algorithms, and hence wavelets are now used in 
both discrete problems, as well as in harmonic analysis decompositions. It is our 
purpose to use tools from systems theory in wavelet problems and also show how 
ideas from wavelet decompositions shed light on factorizations used by engineers. 
Each of the various wavelet families demands a separate class of filters, for the case 
of compactly supported biorthogonal wavelets, see for example Resnikoff, Tian, 
Wells [|60| and Sebert and Zou l63l . By now there is a substantial literature on the 
use of filters in wavelets (see e.g., [18. 26 37 39 1). For filters in wavelets, there are 
two pioneering papers |50]|5l], and the book ll56l . 

In a previous work [;9| we characterized all rational wavelet filters attaining unitary 
values on the unit circle. It turned out that this family is quite small ( and in partic- 
ular the subset of Finite Impulse Response filters, commonly used in engineering). 

Thus, we here remove both restrictions on the filters, i.e., rational and unitary, and 
consider W(z) which are generalized Schur functions, and use reproducing kernel 
Pontryagin spaces associated with W. See J9 for background. 

We hope that this message will be useful to practitioners in their use of these rigor- 
ous mathematics tools. 



3 Pontryagin spaces and Krein spaces 

For a number of problems in the study of signals and filters (for example stability 
considerations), it is necessary to work with Hermitian inner products that are not 
positive definite. This view changes the Hermitian quadratic forms, allowing for 
negative squares, as well as the associated linear spaces. But more importantly, this 
wider setting also necessitates changes in the analysis, for example in the meaning 
of the notion of the adjoint operator, as well as the reproducing kernels. There are a 
number of subtle analytic points involved, as well as a new operator theory. We turn 
to these details below. 
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3.1 Krein spaces 

A Krein space is a pair (V, [•, •]), where V is a linear vector space on C endowed 
with an Hermitian form [•, ■], and with the following properties: V can be written as 
V = V+ + V-, where: 

1. V+ endowed with the Hermitian form [•, •] is a Hilbert space. 

2. V- endowed with the Hermitian form —[•,•] is a Hilbert space. 

3. It holds that V+ n V- = {0}. 

4. For all v± G V±, 

[v+,v_]=0. 

The representation V = V+ + V_ is called a fundamental decomposition, and is 
highly non unique as soon as dim V_ > 0. Given such a decomposition, the map 

(r(v+ + V- ) = v+ — V- 

is called a fundamental symmetry. Note that the space V endowed with the Hermi- 
tian form (where w = w+ + w- is also an element of V, with w± G V±) 

(v, w) = [v, ffw] = [v+,w+] — [V-,W-] 

is a Hilbert space. These norms are called natural norms, and they are all equiva- 
lent. The Hilbert space topologies associated to any two such decompositions are 
equivalent, and V is endowed with any of them; see ITT31 p. 102]. When V_ is finite 
dimensional, V is called a Pontryagin space and the dimension of V- is called the 
negative index (or the index for short) of the Pontryagin space. We refer to the books 
lfl2ll . fl31 . 11361 . iU for more information on Krein and Pontryagin spaces. Note that 
in 11361 it is the space V+ rather than V_ which is assumed finite dimensional in the 
definition of a Pontryagin space. Surveys may be found in for instance in 1301 . 11291 , 
Q- It is interesting to note that Laurent Schwartz introduced independently the no- 
tion of Krein and Pontryagin spaces (he used the terminology Hermitian spaces for 
Krein and Pontryagin spaces) in his paper [62]. For applications of Krein spaces 
to the study of boundary conditions for hyperbolic PDE, including wave equations, 
and exterior domains, see for example [23 52 53 58|. We now give two examples, 
which will be important in the sequel. 

Example 3.1 Let J G C pxp be an Hermitian involution, i.e. 

J = J- 1 =J*. 

Such a matrix is called a signature matrix. We denote by Cy the space C endowed 
with the associated indefinite inner product 

[x,y]j=y* Jx , x,yeC p . 

It is a finite dimensional Pontryagin space. 
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Example 3.2 Let J be a signature matrix. We consider the space H2(B) P of func- 
tions analytic in D and with values in C p : 

f(z) = f j a„z n , a n e£", 

n=0 

such that 

oo 

L U n a n < °°- 
«=0 

Then, H2(B) P endowed with the Hermitian form 

co oo 

\f, 8 ]j = £ M g(z) = £ V") 

is a Krein space, which we denote by H2,./(B). 

In the above example, if p = 1 and J = 1 (as opposed to J = — 1) the space H2,/(B) 
is equal to the classical Hardy space H^B) of the unit disk. 



3.2 Operators in Krein and Pontryagin spaces 

When one considers a bounded operator A between two Krein spaces (Jfr\ ,[•,•] i ) and 
(J%2, [vb) ( m this P a P er > it will be most of the time between two Pontryagin spaces) 
the adjoint can be computed in two different ways, with respect to the Hilbert spaces 
inner products, (and then we use the notation A*) and with respect to the Krein 
spaces inner products (and then we use the notation A M). More precisely, if <7i and 
(72 are fundamental symmetries in and J^2 which define the Hilbert spaces inner 
products 

= and (/2,^2>2 = [02/2,g2]2, 

(with fi,g\ £ and f 2l g2 € J^), we have for f\ £ Jt[ and f 2 6 ,%fi 

[A/l,/ 2 ] 2 = (cT2A/l,/2>2 

= (fi,A*a 2 f 2 h 
= [fiA [ * ] f2]u 

with 

aW=CT!A*C72- (6) 
In the case of Cj (see Example 13. It we have 

AM=/AV. (7) 

The operator A from 2#(A) C where (J^ , [•, -]i) is a Krein space, into the Krein 
space (J^2, [•, ^2) is a contraction if 
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[AkuAhh^lhMh, V*iG^(A). 

A densely defined contraction, or even isometry, operator A between Krein spaces 
need not be continuous, let alone have a continuous extension. See for instance [29, 
Theorem 1 . 1 .7] . In the case of Pontryagin spaces with same negative index, A has a 
continuous extension to all of see |6] Theorem 1.4.1, p. 27], and Theorem [33] 
below. Even when it is continuous and has a well-defined adjoint, this adjoint need 
not be a contraction. The operator is called a bicontraction if both it and its adjoint 
are contractions. When the Krein spaces are Pontryagin spaces with same negative 
index, a contraction is automatically continuous and its adjoint is also a contraction. 
An important notion in the theory of Pontryagin spaces is that of relation. Given 
two Pontryagin spaces 3?\ and a relation is a linear subspace of 3P\ x For 
instance the graph of an operator is a relation. The domain of the relation St is the 
set of f £ such that there is a g £ £?2 for which (f,g) £ Si. K relation 3% is 
called contractive if, 

[g,gh<[f,fh V(/,s)e^. 

A key result is the following theorem of Shmulyan (see |6] Theorem 1.4.1, p. 27]). 

Theorem 3.3 A densely defined contractive relation between Pontryagin spaces 
with same negative index extends to the graph of a uniquely defined contraction 
operator from & > \ into 3?i- 



3.3 Kernels 

Recall that a (say, matrix- valued) function K(z, w) of two variables, defined for z 
and w in a set Q is called a positive definite kernel if it is Hermitian: K(z,w)* = 
K(w,z) for all z, w £ Q, and if for every choice of M £ N and w\, . . . , wm £ Q the 
M x M Hermitian block matrix with (£, j) block entry K(wi,Wj) is non negative. 
For instance, if b is a finite Blaschke product, 



^)=n T 



m z-a n 



n=l 



~a 



for some a\ , . . . , a m in the open unit disk, the kernel 

l-b(z)b(w)* 



h{z,w) 



l-zw* 

is positive definite, as can be seen from the formula 

h{z,w) = (k b (-,w),k h (-,z))n 2 
When b is replaced with a function s analytic and contractive in the open unit disk, 

l-zw* 



the corresponding kernel k s (z, w) = 1 is still positive definite in D, see |[T6l . 
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IfrTll . This follows, for instance, from the fact that the operator of multiplication by 
s is a contraction from H2(B) into itself. In the special case of a finite Blasckhe 
product (or more generally, of an inner function), this multiplication operator is an 
isometry. This makes the underlying computations much easier. More generally, the 
kernels which appear in the following section can be seen as far reaching general- 
izations of the kernels kb(z,w). 

The notion of positive definite kernel has been extended by Krein as follows: 

Definition 3.4 Let k e No- A ( say, matrix-valued) function K(z, w) defined on a 
set Q has K negative squares if it is Hermitian, and if for every choice of M E 
N and w\, . . . ,wm E Q the M x M Hermitian block matrix with block entry 

K(wi,Wj) has at most K strictly negative eigenvalues, and exactly K strictly negative 
eigenvalues for some choice ofM, w\, . . . , wm- When K = 0, the function is positive 
definite. 

The one-to-one correspondence between positive definite kernels and reproducing 
kernel Hilbert spaces was first extended to the indefinite case by L. Schwartz; see 
l62l : There is a one-to-one correspondence between reproducing kernel Pontryagin 
spaces and kernels with a finite number of negative squares. For completeness, we 
mention that such a result fails if the number of negative squares is not finite. A 
necessary and sufficient condition for a function to be the reproducing kernel of a 
Krein space is that this function is the difference of two positive kernels, but the 
associated Krein space need not be unique. Here too we refer to Schwartz l62l . and 
also to the paper [1]. Realization of operator- valued analytic functions (without as- 
sumptions on an associated kernel, but with some symmetry hypothesis) have also 
been considered. See for instance 11271 . The C pxp -valued function K(z,w) defined 
for z,w in an open set £2 of the complex plane will be called an analytic kernel if 
it is Hermitian and if it is analytic in z and w* . If it has moreover a finite number 
of negative squares, the elements of the associated reproducing kernel Pontryagin 
space are analytic in Q.. See |6| Theorem 1.1.2, p. 7]. 

There are two important classes of operators between reproducing kernel spaces, 
namely multiplication and composition operators. We conclude this section with 
three results on these operators. 

Theorem 3.5 Let (J^i ,[•,•] l ) and (^2,[vk) be two reproducing kernel Krein 
spaces of vector-valued functions, defined in Q, and with reproducing kernels 
Ki(z,w) andK 2 (z,w), respectively C PlXpi - and C P2Xpi -valued. Let m be a C 2 ^ 1 - 
valued function and let <p be a map from Q into itself. Assume that the map 

(T m , 9 f)(z)=m(z)f((p(z)) (8) 

defines a bounded operator from (J^i,[v]i) into (J^2, [•, -]2)- Then, for every z,w € 
Q, and & G C P2 , 

(T^K 2 (;w)^ (z)=K l (z,( P (w))m(w)*l; 2 . (9) 
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Proof: Let z, w G Q, t, 2 G and ^ G C p ' . We have 

li (rit^ 2 (-,w)&) (z) = [rW^aO.wj&.ZiO.z^ili 
= [z 2 (.,w)& J r m , 9 (z 1 (.,z)§ 1 )] 2 

= [K 2 (;w)&,m(-)K 1 ((p(-),z)Zi}2 
= [m(-)K l (( P (-),z)$i,K 2 (-,w)&]* 2 

= (^m(w)K 1 ( ( p(w),z)^T 
= § 1 *Zi(z,9(w))m(w)*§ l . 



As a corollary we have the following result: 

Theorem 3.6 Assume in the preceding theorem that and ,J€ 2 are Pontryagin 
spaces with same negative index. Then, T m m is a contraction if and only if the kernel 

K 2 (z, w) - m(z)^i (p(w))m{w)* (10) 

is positive definite in Q. 

Proof: Assume that T is a contraction. Then, its adjoint is also a contraction since 
the Pontryagin spaces have the same negative index. Let g G be of the form 

N 

g(z) = K 2 (z,w k )£ k , 

k=l 

where N G N, W\ , . . . , wjy G CO and §i , . . . , ^ N G C' 2 . By © we have 

N 

t,k=i 

N N 

= [£ K t (z, <p(w k ))m(w k )*% k , £ K x (z, 9(w/))m(w/)*§<] i 

k=l t = \ 

_ \rp[*] rp[*] 1 

— l-'nijip^i lm,<p8\ 1 

2V 

= XI &K 2 (we,w k )% k , 

£,k=l 

and hence the kernel dTOb is positive definite. Conversely, assume that the kernel 
dTOb is positive definite. Then the linear span of the pairs of functions 

(tf 2 (.,v^,^(-,(pW)mW^), w G £2,1; G C« 
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defines a linear densely defined contractive relation in J(f\ x Jfa, By Shmulyan's 
theorem (see Theorem |3.3t , this relation has an everywhere defined extension which 
is the graph of a bounded contraction: There is a unique contraction X from 3fi into 
such that 

X(K 2 (-,w)%)=K 1 (;(p(w))m(w)*Z, we £2, £,eC P2 . 
By (O, we have X W = T m ^, and this concludes the proof. □ 

We will consider in the sequel special cases of this result, in particular when 

m(z) = (1 z ■ ■ ■ z N ~ l ) , 
see Theorem l5.3l or more generally when 

m(z) = (m (z) m\ (z) ■■■ m N _ x (z)) , 

see Theorem 15.41 The operator T m ,m defined by ([8]l is then a block operator, and 
its components satisfy, under appropriate supplementary hypothesis, the Cuntz re- 
lations formally defined in (fl~8T> - (fT~9l> below. 

We conclude this section with a result on composition operators in reproducing 
kernel Pontryagin spaces. 

Theorem 3.7 Let K(z,w) be a <C px p -valued function which has K negative squares 
in the set Q. The associated reproducing kernel Pontryagin space will be denoted 
by &[K). Let <p be a map from Q into itself and assume that: 

/(p(z)) = =*>/=<) 

forfe &>{K). Then: 

(a) The function Kq,(z,w) = K((p(z),(p(w)) has at most K negative squares in Q 
and its associated reproducing Pontryagin space is the set of functions of the form 
F(z) = f((p(z)), with f G tP(K) and Hermitian form 

[ F ,G]^(K 9 ) = \f,g]^{K)- (11) 

(2) The map f i— > /(<p) is unitary from ZP{K) into itself if and only if 

K(z,w)=K(q>(z),q>(w)), Vz,we£2. (12) 

Proof: Set 

.* v = {f(9(z)),fe.nK)}- 

By hypothesis, we have /(<p(z)) = if and only if f = Q, and so the Hermitian form 
(fTTT i is well defined and induces a Pontryagin structure on Furthermore, with 

c G C p and F(z) = f(<p(z)) G JK<p, we have: 
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[f(-),^(-,w) c ]^ w = [/(•),*(•, <pM>W) 

= c*/(9W) 
= F(w>), 

and hence the reproducing kernel property is in force. To prove (b) we use the 
uniqueness of the kernel for a given reproducing kernel Pontryagin space. □ 

To fine-tune the previous result, note that for (p(z) = z N , the composition map is 
an isometry from H2(H>) into itself, but is not unitary (unless N = 1). We also note 
that the preceding theorem holds also for reproducing kernel Krein spaces. Indeed, 
the correspondence between functions which are difference of positive functions 
on a given set and reproducing kernel Krein spaces is not one-to-one, but a given 
reproducing kernel Krein space has a unique reproducing kernel. 



4 Generalized Schur functions and associated spaces 

In this section we review the main aspects of the the realization theory of generalized 
Schur functions and of their associated reproducing kernel Pontryagin spaces. 



4.1 Generalized Schur functions 

In the positive definite case, this theory originates with the works of de Branges 
and Rovnyak, see |[T6l[T7l . In earlier work on models involving operators in Hilbert 
space, and matrix factorization, de Branges spaces have served as a surprisingly 
powerful tool. The theory was developed in the indefinite case in in a fundamental 
series of papers by Krein and Langer, see for instance ll44l l45l l46l l47l l48l . and, 
using reproducing kernel methods in JTJ and in the book |6 j. It was later used in 
(6] p. 119] and in the paper J3] to study generalized Schur functions with some 
given symmetry. In this paper we use this setting to present non rational and non 
unitary wavelet filters. In (6) the case of operator valued functions is studied, but 
we here consider the case of C pxp -valued functions. We now recall the definition 
of a generalized Schur function. A (say C /,x/, -valued) function W is called a Schur 
function if it is analytic and contractive in the open unit disk, or, equivalently, if the 
associated kernel 

I p -W(z)W(w)* 
Kw(z,w) = — ; (13) 

is positive definite in a neighborhood of the origin. Then, it has a unique analytic 
extension to the open unit disk, and this extension is such that the kernel Kw is 
still positive definite in ED. There are two other kernels associated to W, namely the 

kernel %(z,w) (with W(z) = W(z*)*), and the kernel 
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D w (z,w) 



Wlz)-W(w*) „ i \ 



z— w* 



These three kernels are simultaneously positive definite in the open unit disk. The 
first is the state space for a unique coisometric realization of W, the second is the 
state space for a unique isometric realization of W, and the reproducing kernel 
Hilbert space with reproducing kernel Dw is the state space for a unique unitary 
realization of W. In these three cases, uniqueness is up to an invertible similarity 
operator. 

Let J £ C pxp be a signature matrix. We now consider functions with values in Cj 
defined in Example 13. II denoted by (rather than W). A C /7X/7 -valued functions 
analytic in a neighborhood of the origin is called /-contractive if the associated 
kernel 

K j-m^y 

1 — ZW* 

is positive definite. It has a unique meromorphic extension to the open unit disk, 
and this extension is such that the kernel K@ is still positive definite in the domain 
of analyticity of in D. Here too, besides the kernel K@ we have the kernel (z, w) 
and the kernel 

J0(z)-J@(w*)\ 



D {z,w) 



K @ {z z w) 

0(z)J-@(w*)J 



We note that the kernel K@ can be written as 



I p -0{z)Q(wf 
K @ {z,w) = - 



- zw 

where [*] denotes the adjoint in Cj. This conforms with the way these kernels and 
the two other related kernels are written down in J6). 

As we already mentioned, Krein and Langer developed in B4ll45l l46 47 48], the 
theory of operator-valued functions such that the corresponding kernels K@ (with 
a signature operator rather than a signature matrix) has a finite number of negative 
squares in some open subset of the open unit disk. Then, has a unique mero- 
morphic extension to the open unit disk, and this extension is such that K@ has the 
same number of negative squares in Q(&), the domain of analyticity in D. The 
three kernels have simultaneously the same number of negative squares, and as in 
the positive definite case, are respectively state spaces for coisometric, isometric and 
unitary realizations of 0. 

In the special case 7 = 7 (we return to the notation W rather than for the function), 
Krein and Langer proved, see l44l . that W can be written as WqBq 1 , where Wo is 
analytic and contractive in the open unit disk, and where Bq is a finite matrix-valued 
Blaschke product. It follows that W has a finite number of poles in the open unit 
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disk. In the rational case, and when W takes unitary values on the unit circle, W is 
a quotient of two matrix-valued rational Blaschke product. Note however that when 
/ has mixed inertia, W may have an infinite number of poles, even when k = 0. For 
example, take 

/= (J-°i) and *«=G*«-0' 

where b is a convergent Blasckhe product with an infinite number of zeros. Such 
examples originate with the work of Potapov |59l . 

Definition 4.1 We denote by S*g xp (p) the family of C pxp -valued functions W 
meromorphic in the open unit disk, and such that the kernel Kyy (defined by ( 113l >) 
has K negative squares in the domain of analyticity ofW in D. 
Given a signature matrix J, we denote by S*£(J$) the family of C pxp -valued func- 
tions © meromorphic in the open unit disk, and such that the kernel K@ (defined by 
d 14b J has K negative squares in the domain of analyticity of & in D. 
Wedenoteby £P(W) and 3?(&) respectively the associated reproducing kernel Pon- 
tryagin spaces. 

Since the kernels Kw and K® are analytic in z and w*, the elements of the associated 
reproducing kernel Pontryagin spaces are analytic in the domain of definition of W 
or respectively. See (6, Theorem 1.1.3, p. 7]. 

More generally, it is useful to consider non square generalized Schur functions. We 
consider J\ € C p,xp ' and J2 £ C Pl xpi two signature matrices, of possibly different 
sizes, such that <(3j is in form denoted by v_ : 

V_(Ji) = V_(J 2 ). 

Reproducing kernel Pontryagin spaces with reproducing kernel of the form ©: 

J 2 -&{z)Ji@{w)* 
1 ~zw* 

when @ is C P2Xpi -valued and analytic in a neighborhood of the origin, have been 
characterized in (6J Theorem 3.1.2, p. 85] (in fact, the result there is more general 
and considers operator- valued functions). In the statement below Rq denotes the 
backward-shift operator 

z 

Theorem 4.2 Let [■,■]&) be a reproducing kernel Pontryagin space of (De- 
valued functions. It has a reproducing kernel of the form © and only if it is 
invariant under the backward-shift operator Rq and 

k<kI r < [fJ}&-Aoyj 2 f(o), v/ g 9. 
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4.2 State spaces and realizations 

We begin with recalling the following definition. Let W be an operator- valued func- 
tion analytic in a neighborhood of the origin. A realization of W is an expression of 
the form 



where D = W(0) and A,B,C are operators between appropriate spaces. It is an im- 
portant problem to connect the properties of W and of the operator matrix 



When the values of W are linear bounded operators between two Krein spaces, Az- 
izov proved that a realization exists, and that M can be chosen unitary. See ifTTI . 
and see l27l for further discussion and additional references. When W is a matrix- 
valued rational function without a pole at the origin, the spaces may be chosen finite 
dimensional, when no special structure is forced on the operator matrix M. 

In Section 14.11 we have studied the correspondence between kernels and operator 
valued Schur functions. Here we then pass to the realizations of Schur functions. The 
introduction of Schur functions offers many advantages, relevant to algorithms and 
to computation. Case in point: In the next subsection, we give explicit formulas for 
realizations, i.e., for the computation of the four block operator entries A through 
D making up admissible realizations of a given Schur function, and therefore of 
a kernel. As we show, there are several such choices, the coisometric realization 
(Theorem 14. 31 l. and the unitary realization of de Branges and Rovnyak (Theorem 
14.4b , among others. There is in turn a rich literature on Schur algorithms in various 
special cases, see for example Q for an overview. In preparation of Section l4~3l we 
need some definitions. Let 3? denote the space where A acts in ( TT5T >. We say that the 
realization is closely inner connected if the span of the functions 



where E, runs through C p (recall that J £ C pxp ) and z runs through a neighborhood 
of the origin, is dense in With the same choices of t, and z, it will be called 
closely outer connected if the span of the functions 



W(z) =D + zC(I-zA)- 1 B, 



(15) 




(16) 



is dense in J 2 , and connected if the span of the functions 
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(I-zA)- l B%, and (/- wA^^C^T] 

is dense in 3? (77 running through CP and w through the same neighborhood of the 
origin as z). Here the adjoints are between Pontryagin spaces. We note that the ter- 
minology is different from that of classical system theory. In the finite dimensional 
case, what is called here closely inner connected corresponds to observability, and 
what is called outer connected corresponds to controlabilty. The notion of being 
closely connected is specific to this domain, and is, in general, different from mini- 
mality. 



4.3 Coisometric and unitary realizations 

Let © G be a generalized Schur function, assumed analytic in a neighborhood 
of the origin. In this section we review how the spaces 3P(®~) and 2>(&) are the 
state spaces for coisometric and unitary realizations respectively. For the following 
theorems, see J6j Theorem 2.2.1, p. 49] and Theorem 2.1.3] respectively. In 
Theorems 14. 3l and 14. 4l below the notions of coisometry and unitarity means that M 
in ( fT6l ) is an operator coisometric (resp. unitary) from the Pontryagin ^(@)©C/ 
into itself (resp. from S>(0) © Cj into itself). 

Theorem 4.3 Let J G C pxp be a signature matrix, and & € S?^ be analytic in a 
neighborhood of the origin. Then the formulas 

A/(z) = /WzM, 

z 

,„. w , ©(z)-©(0)„ 



Cf = f(0), 
D$=0(O)£, 

with f G &(&) and E, G C, define a closely outer connected realization of@ which 
is coisometric. This realization is unique up to a continuous and continuously invert- 
ible similarity operator. 

This coisometric realization was introduced by L. de Branges and J. Rovnyak in 
lfl6ll for scalar Schur functions, and extended to the operator- valued case in ifTTll . 
We note that the coisometric realization is also known as the backward shift realiza- 
tion; see e.g. ll33ll . 

L. de Branges and J. Rovnyak also formulated the unitary realization below. 

Theorem 4.4 Let J G C pxp be a signature matrix, and G J?^ be analytic in a 
neighborhood of the origin. The formulas 
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zg(z)-0(z)Jf(O), 

e{z)-e(o) 

(J-0(z)J&(oy)£ 



c ({)=/(°)> 



D§=©(0)§, 

with f € and E, € C p , define a closely connected realization of which is 

unitary. This realization is unique up to a continuous and continuously invertible 
similarity operator. 

It is important to note that, in some cases, all three realizations are unitary. This is in 
particular the case when is rational and /-unitary on the unit circle. See Section 
S3] 



4.4 Finite dimensional de Branges spaces 

The finite dimensional case is of special importance, and the case J = I was con- 
sidered in details in our previous work (9). Then the three realizations are unitary, 
and it is easier to focus on the &(0) spaces. As proved in [7, Corollary p. 1 1 1] for 
the case J = I and in j7. Theorem 5.5, p. 1 12] for the general case, given € J^, 
the associated space £P(0) is finite dimensional if and only if is rational and J 
unitary on the unit circle: 

0{e")*J0{e")=J, 

at all points e" (t 6 [0,27r]) where it is defined. If moreover is analytic in a neigh- 
borhood of the closed unit disk, we have 

^(0)=H 2J e0H 2J . 

Rationality is not enough to insure that ,<3?(0) is finite dimensional, as illustrated 
by the case J = 1 and = 0. Then, &>{0) = H 2 (D). 



Definition 4.5 We will denote by the multiplicative group of rational C pxp - 
valued functions which take J -unitary values on the unit circle, and for which the 
corresponding kernel K@ has K negative squares. We set 

w J = U ^i. 

K=0 
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The results and realizations presented in the previous section take now an easier 
form. The various operators can be seen as matrices. Unitarity above is with respect 
to the indefinite metric of ^(©)(BCj, and we can rephrase Theorem l4~4l as follows: 

Theorem 4.6 Let W be a rational <C px p -valued function analytic at the origin, and 
let 

W(z)=D + zC(I-zA)- l B 

be a minimal realization ofW. Then, W is J -unitary on the unit circle if and only 
if there exists an invertible Hermitian matrix H (which is uniquely determined from 
the given realization) such that 

(^yn(c£Ho°> 

The change of variable zh> 1/z yields: 

Theorem 4.7 Let W be analytic at infinity, and let 

W{z)=D + C{zI-A)- l B. 

be a minimal realization of W. Then, W is J -unitary on the unit circle if and only 
if there exists an invertible Hermitian matrix H (which is uniquely determined from 
the given realization ) and such that (1 1 "7b holds. 

The matrix H is called the associated Hermitian matrix (to the given minimal real- 
ization). This result was proved in (SJ Theorem 3.10] for the case where A is non- 
singular. For the approach using reproducing kernel Hilbert spaces, see EH 5]. 



5 Cuntz relations 



5.1 Cuntz relations and the de Branges-Rovnyak spaces 

The results of this section are related to l22l and iflOll . In that last paper, the func- 
tions 1, . . . ,z N ~ l below are replaced by the span of a finite dimensional backward- 
shift invariant subspace, but the discussion is restricted to the Hilbert space case and 
scalar-valued functions. 

Normally by Cuntz relations we refer to a finite system of isometries Si , . . . , Sn in a 
Hilbert space Jf satisfying two conditions: 

(a) Different isometries in the system must have orthogonal ranges, 



S*jS k = 0, ;V k, 



(18) 
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and 

(b) The sum of the ranges equals 

N 

ZSjS*=I^. (19) 

.7=1 

Note that (a) already forces to be infinite dimensional. Indeed, if Jtff is finite 
dimensional, an isometry is unitary and the orthogonality of the ranges is not possi- 
ble, see the discussion below and Section [531 If we allow the isometries to operate 
between two finite dimensional spaces of different dimensions, then one can find 
isometries which satisfy the Cuntz relations. It is the set of three conditions: Each 
Si is isometric in a Hilbert space J^, and (a) and (b), together imply that every 
realization is a representation of a simple, purely infinite C*-algebra, called On- In 
applications to filters, the N individual subspaces represent frequency bands. This 
allows for versatile computational algorithms tailored to multiscale problems such 
as wavelet decompositions, and analysis on fractals. In our present paper, we relax 
some of the original very restrictive axioms, while maintaining the computational 
favorable properties. Our more general framework still allows for algorithms based 
on iteration of the operator family Si, . . . ,Sft in a particular representation. 

If one allows isometries between two Hilbert spaces, then the finite dimensional 
case may occur, as illustrated by the following example: 

J^i = C, J? 2 = C 2 , 




We have 

S\Si = S* 2 S 2 = 1 , S\S 2 = S* 2 Si = 0, 

and 



SiS\+S 2 S* 2 = 




We go beyond the setting of Hilbert space, and relax the conditions (a) and (b) 
imposed in the original framework from C* -algebra theory, allowing here isomet- 
ric operators between two Pontryagin spaces. We still preserve the features of the 
representations of use in iterative algorithms. 

It is not surprising that in Section 1531 we have finite dimensional spaces. Now for 
the generalized theory, we must allow for de Branges and Rovnyak spaces, and for 
negative squares and signature matrix. The resulting modifications in the form of 
the Cuntz relations, in the case of Hilbert space, entails some non-trivial modifica- 
tions addressed in the next two sections. Our main results for this are proved in the 
present section, and in Section l531 for the finite dimensional case. 
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The main result of this section is that one can associate in a natural way to an element 
£ ,5^(D) a family of operators which satisfy the Cuntz relations. We begin with 
a preliminary result, which is a corollary of Theorem l3.7l with (p(z) = z N ■ 

Proposition 5.1 Let & <G .5^(0), and let be the associated Pontryagin 

space, with reproducing kernel 

K & (z, w) = 



1 — zw* 
The function 

„ „ j-&{z N )j©{ W N r 

W ,W ) = X _ Z N W *N 

has also K negative squares in its domain of definition in D. The associated repro- 
ducing kernel Pontryagin space „#jy is equal to the space of functions of the form 
F(z) = f(z N ), where f G with the following indefinite inner product 

[F,G]^ N = \f,g]&(0), (20) 

where g G &>{&) and G(z) = g(z N ). 
We have: 

Theorem 5.2 Let & G J?^(1D)), and let be the associated Pontryagin space 

with reproducing kernel 

A@ [Z, W) 



I — zw* 

Then, for N E N, the function ©n defined by &n(z) = &(z N ) belongs to S^i K . Fur- 
thermore, consists of all the functions of the form 

Any such representation is unique, and the inner product in is given by 

N-l 

\f,g]0»(e N ) = 2- \fhS ]]&>{&)■> 

7=0 

where g{z) = Y!J=o z J gj(z N )for some go,... ,gjv_i G 
Proof: We proceed in a number of steps. 

STEP 1 : It holds that V-{0 N )<N- K. 



Indeed, 
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j-e(f)j0^y _ j-0(z n )j&{w n )* i-zV* 

l-zw* ~ l-z N w N * l-zw* 



J-0(z N )J0(w N ) * 
l-z¥' 



(E z 

k=0 



This expresses the kernel K@ N as the sum of N kernels, each with K negative squares. 
Thus, V-(&n) < Nk. To show that there is equality, we need to show that the asso- 
ciated spaces have pairwise intersections which all reduce to the zero function. 

STEP 2: Letk,£e {0,...,N- 1}, such thatk^ I. Then, with Jt N as in theprevious 
theorem: 

z k ~£ N r)z ( ^ N = {0}. 

Indeed, assume that k > £ and let f,g <G be such that 

z f{z )=Zg{Z ). 

Then, / and g will simultaneously be identically equal to Q px \. Assume / ^ px i- 
One of its components, say the first, with / = (*i(z) • • • Xp(z))' is not identically 
equal to zero (p is the size of the signature matrix J). Then we obtain 

xi(z N Y 

where y\ denotes the first component of g. Since / and g are meromorphic in D, 
the function yi/x\ has a Laurent expansion at the origin. Moreover the Laurent 
expansion of ^(z N ) contains only powers which are multiple of N. By the unique- 
ness of the Laurent expansion, this contradicts the fact that it is equal to z k ~ t , with 

\k-£\<N. 

STEP 3: It holds that 

and it holds that v@ N = Nk. 

This is because the spaces z } .^n have pairwise intersections which reduce to the 
zero functions in view of STEP 2. □ 



Theorem 5.3 In the notation above, set 

(Sjf)(z)=z J f(z N ) &{0)^&{0 N ). 

Then, 

Sff = fj &(0n)^&(@), (21) 

and 
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^ [*] _ ( 22 ) 

7=0 

where the [*] denotes adjoint between Pontryagin spaces. 
Proof: We proceed in a number of steps. 

STEP 1: The operators Sj are continuous. 

The operators Sj are between Pontryagin spaces of different indices, and some care 

is required to check continuity. To this end, fix j G {0, ,iV — 1} and note that Sj 

is everywhere defined. Furthermore we claim that it is a closed operator. Indeed, 
let fi ,/2 . . . be a sequence of elements in &(&) converging strongly to / G IP(0) 
and such that the sequence Sjfi , Sj/2 , . . . converges strongly to g G 3? (0;v) . Strong 
convergence in a Pontryagin space implies weak convergence, and in a reproducing 
kernel Pontryagin space, weak convergence implies pointwise convergence. There- 
fore, for every w where is defined, 



and 



lim/*(w) =/(w), 



lim (Sjf k )(w)=g(w). 

k— >°o 



Since (Sjfk){w) = w J fk{w), and thus g(w) = w j f(w). Therefore g = Sjf, and the 
operator Sj is closed, and hence continuous. 

STEP 2: (EB is in force. 

Letg(z) =Zk=o z k 8k{z N ) G 0>{©n) where the^ G ^>(0), and let u G ^(0). Then, 

AT— 1 
k=0 

= [u, Sj g]s»{&)> 

where [, ]&>(&) and [, ]g>(@ N ) denote the indefinite inner products in the correspond- 

[*] 

ing spaces. Hence, we have S) g = gj. 
STEP 3: The Cuntz relations hold. 



From (|2TT > we have for u £ &(©) 
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Furthermore, for f(z) = YJJ^q z j fj(z N ) € ^(©n) (where the fj G &>(®)\ we have 



if j?k, 
u if j = k. 



and thus 



s k st ] f = s k (f k ) = z k f k (z N ), 



N ~ l [*] 

A-=0 



□ 

We note that, with 

S = (So Si — S N -i) @>{&) N — >• 
the Cuntz relations d2TT l can be rewritten as 

SSM=/^ (0jy) and sMs = /^ (0) *. 
At this stage, let us introduce some more notation. We set 

® Nk {z) = &{z Nk l 

and sj ' = Si for i = 0, . . . ,N — 1. We can reiterate the preceding analysis with 0^ 

instead of 0. We then obtain AT isometries S^\. . . (Siilj from ^(0^) into &(© N 2) 
satisfying the Cuntz relations. Iterating A: times, one obtains k sets of isometries, 

oCM) cO'-l) fc 

from &{@ N j-\) into £P(© N j), which also satisfy the Cuntz relations. This gives us 
A^' isometries 

o(0)eW...e(n) 

h h h ' 

with (h,h,---,h) G {0,...,AT- 1}*, from ^(0) into ^(0 N k), all satisfying the 
Cuntz relations. 



5.2 Cuntz relation: The general case 

We now wish to extend the results of Section 15.11 and in particular Theorem 
to the case where the N functions l,z, . . . ,zr are replaced by prescribed func- 
tions mo(z),mi (z), . . . ,Wjv-i(z), whose finite dimensional linear span we denote by 
and the kernel K@{z,w) is replaced by a given analytic C A ' xA '-valued kernel 



24 



Daniel Alpay, Palle Jorgensen and Izchak Lewkowicz 



K(z,w) and the kernel K@ N {z,w) is replaced by a kernel K(z,w). Let as in Section 
15.11 Kn(z,w) = K(z N ,w"). We address the following problem: Given K and K two 
Hermitian kernels defined on a set £2, and with a finite number of negative squares 
there, when can one find decompositions of the form 



f(z)= N fm n (z)g n (z" 

n=0 



(23) 



where the functions go, . . . ,g^-i belong to Z?(K) for some, or all, elements in 
&>(K). We have: 

Theorem 5.4 Let K(z, w) and K(z, w) be two kernels defined on a set Q, and as- 
sume that 

V-{K)=NV-(K). (24) 
Let mo, . . . ,ffJAT-i be N functions on Q. Assume that the kernel 

_ N-l 

K{z,w) - ( £ m n (z)m n (w)*)K{z,w) 

is positive definite in Q. Then, with <p{z) = z N , the choice g„ = T^^f n , n = 
0,...,N—1 solves I 



Proof: We use Theorem l3.6l with K2(z,w) = K(z,w) and 



Ki(z,w) 



(K{z,w) ••• \ 
K(z,w) ••• 



K(z,w)J 



Then 



N-l 



( Y, mn(z)m n (w)*)K(z,w) = m(z)K\ (z)m{w)* , 

n=Q 

and Theorem |3.6| with K\ and K-i as above, and 

m(z) = (m (z) mi (z) ■■■ m#-i (z)) , and <p(z) = z N , 
leads to the fact that the map 

f^m(z)f(z N ) 

is a contraction from (&>(K)) N into 2? (K ) . □ 

In applications, one uses the kernel K(z,w) = Kn(z,w) in the above result. 
Proposition 5.5 A sufficient condition for ( 124b to hold is that 



mj^>(K) n m k ^(K) = {0} , 



(25) 
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for all j,k G {0, . . . ,N — 1} such that j ^ k. 

Proof: Indeed, when this condition is in force, we have that the Pontryagin space 
with reproducing kernel m(z)K\(z,w)m(w)* is the direct sum of the Pontryagin 
spaces with reproducing kernels mj(z)K(z,w)mj(w)* , j — 0, . . . ,N — 1. □ 

We note that there are similarity between d23l l and the solution of Gleason's prob- 
lem: Gleason's problem is the following: Given a linear space of functions ^ of 
functions analytic in a set Q C C , and given a € £2, Gleason's problem is the fol- 
lowing: when can we find functions gi,---,gN £ ^# (which depend on a) such that 

N 

f(z)-f{a) = Y,{z n -a n )g n (z,a) 

n=\ 



5.3 Cuntz relations: Realizations in the rational case 

Recall that for a given generalized Schur function 0, we presented in Theorems 
l4.3l and l4~4l coisometric and unitary realizations respectively. The unitary realization 
turns to be more involved than the coisometric backwards shift realization. In some 
cases, these two realizations are unitarily equivalent, in particular when is ratio- 
nal and /-unitary on the unit circle. As we already discussed in Section H~4l this is 
equivalent to having the space £?(&) finite dimensional. In this section we adopt 
this simplifying assumption and study the realization of ©n(z) = & (z N ) in terms of 
the realization of ©. 

We take the signature matrix / to belong to C LxL . We know (see 0[6] and Theorem 
@3] above) that 

&(z N ) = ^ + z c S'{I-z^y l 3§ 

where & = ®/v(0) = 0(0) and where and c £ are defined as follows: c & is the 

evaluation at the origin, 

«y=/(o). 

3& is defined by 

%( Z )-0„(O) ^ 

z 

and srf is the backward shift in ^(0^). The matrix, see 0, 



is unitary in the 3^(0^) metric. We know from Theorem 15. 21 that is equal 

to the space of functions of the form 
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N-i 



/(z)=£zV*(A 
where the E &(©) are uniquely defined. We will denote by U the map 



(26) 



ffo\ 
h 



from £P(©n) onto (&(@)) N . In view of d22l) . U is a unitary map (between Pontrya- 
gin spaces). 

Let T denote the following map from (&{&)) N into itself defined by 



TUf: 



/0 / ■ • • 0\ 
/ ■•• 



\RoOO--Oj 



Vn-iJ 



h 



VW 



Proposition 5.6 Let f e &(©n), with representation d26l l. It holds that 

U £ /f=(TUf)(z N ) 1 (27) 

and it holds that 

(^f^g) n@N) = (TUfJUg) {;n0))N . (28) 
Proof: Indeed, with / is of the form d26| i, we have 



k=\ 



so that Ug/U*f is equal to 



/ /o \ 






/l 




h 















that is ( f27T > in in force. Finally d28l ) follows from the formula for the inner product 
in &>{® N ). □ 
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Proposition 5.7 Let f e s {©n) with representation d26l ). Then, 

?/ = C(7 z 0-0) Uf, 
where C is the evaluation at the origin in 3 s (®). 
Proof: This is clear from 



tf/ = /o(0) = (COO ■■■<)) 



h 



Vn-J 



a 

Proposition 5.8 We have 

.^=z N - x {B^){z N ) 
where B is the operator from C L into &(&): 

and we have 

(aS,aTi)p m = (B^Br,)^ {0) , r],$ g C L . 

Proof: We have 

0(^-0(0) 



?$(z)=R o N t;(z) = 



27 



(29) 



(30) 



Equality d30l ) follows form the definition of the inner product in ^(0^). □ 

These various formulas allow to show directly that the realization is indeed unitary, 
and to compute the associated Hermitian matrix in the finite dimensional case. 



6 Decompositions 

6.1 Generalized down-sampling and an Hermitian form 

In the preceding section we considered decompositions of a function in the form 
( |23T >. Here we consider different kind of decompositions. We consider matrices P G 
C NxN satisfying 
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dst(I N -e e N P e )^0, £ = l,...,N-l, and P N = I N - (31) 

We do not assume that P N ~ 1 ^ In, and in particular the choice P = In is allowed. 
The special case P — EnPn plays also an important role. 

Theorem 6.1 Let W be a <C NxM -valued function defined in the open unit disk ( typi- 
cally, M=lorM=N).LetPe C NxN satisfying <ED> and let, for k = 0,. . . ,N - 1, 

Wk(z) = ^ f N f(e N P) ke W(e N z). (32) 



Then, 



k=0 

Proof: We have 

N-l 



W k {e NZ ) = (s N P)- k (W k (z)), k = 0,...,N-l, (33) 
W(z) = £ W k (z). (34) 



W k (e N z) = - £ {e N P) u W{e N e N z) 



N e=0 



(^P)- k (^ N f(^ p ) k{M)w (^^) 



= (e N P)- k (W k (z)), 
since (ENP) kN = In, and this proves (l33l l. To prove (l34b we write 



JV-l JV-l / i N-l \ 

jfe=0 ife=0 \ v £=0 / 

, /iV-1 (N-l \ 

v V=0 \£=0 / 



= W(z), 

since, in view of (l3~Tl i. 

at- 1 

E 

□ 



E = 



AT, if I = 0, 

(/* - (e N P) m )(I - (e N Py y l = Oift = 1 , 2, . . .N - 1 . 



When P = In, the index k = 1 corresponds to the down-sampling operator. 
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In some cases the decomposition d34l l is orthogonal for the underlying Krein space 
(or Pontryagin space) structure. We will assume that the Krein space (J?T, [•, \x) 
consists of C N -valued functions and satisfies the following property: 

Hypothesis 6.2 Let P be a matrix satisfying ( Bil l, and let (p(z) = £nz. We assume 
that: 

(1) The composition operator f t— > f((p) is continuous and unitary from into 
itself. 

(2) The operator of multiplication by P on the left is continuous and unitary from 
Jlf into itself. 

We note that, in particular, the operator Tp^ defined by ©, 

Tp, V f{z)=Pf{ENZ), 
is continuous and unitary from ,J€ into itself. Note also that 

T p,(t> = 

Hypothesis 16.21 hold in particular for the spaces H2,/ when P is /-unitary, that is, 
satisfies 

P*JP = J. 



Theorem 6.3 Let (Jf, [•, •]) be a Krein space of C N -valued functions, satisfying 
Hvpothesis \6.2\ Let W £ J€ and let 

Wk(z) = ^ N t(e N P) ki W(ef f z). (35) 

Then, 

[w h w k ] = o, i + k 

W(z)=W (z) + ---+W N - 1 (z), 

and 

W k (e N z) = (£ N Py k W(z). 

Proof: The last two claims are proved in Theorem 16. II The first claim takes into 
account the hypothesis on J(f, and is proved as follows: We take k{ and ki in 

{0, ,JV — 1}, and assume that £2 <k\. Taking into account the definition of Wk, 

we see that the inner product [W^ , Wj 2 ] x is a sum °f N 2 inner products, namely 

[(e N P) k ^W(e e N l z),(e N P) k ^W(e%z)}jr, t u h e {0,...,N-1}. 

These A^ 2 inner products can be rearranged as N sums of inner product, each sum 
being equal to 0. Indeed, consider first the inner products corresponding to £\ = In. 
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In view of the unitary of the operator Tp q we have 

N f[(e N P) k ^W(e^zUe N P) k ^W{e^z)U = [f ^(e^-pA W,W]x 

= 0. 

Indeed, using < k\ — ki < N— 1, and so, by hypothesis on P, we have 

det(I N -(£ N P) k ^ h -)^l, 

and the sum 

f,=0 

Let us now regroup the factors of \W (z) ,W (£nz)]x ■ Taking into account that 

[P^-i) W (£%- l z),W(z)U = {P kliN - l) W(z),W(e N z)U, 

we have 

N ^[{e N P) k ^W{eiz)M N P) kl - [M) W{£ l N + \)] x + 



+ {(e N P)^ N -Vw(e%- l z),W(z)U 
: [^(e N pf^ e+l ^ + (e N P)^ N - l) ^ W(z),W(e N z)U 

- [{e N P)- h (^{{e N P) h - h )^j W(z),W(e N z)W 
0. 

The remaining terms are summed up to in the same way. □ 



6.3 Decompositions in reproducing kernel spaces 

We begin with a result in the setting of Schur functions, as opposed to generalized 
Schur functions. 

Theorem 6.4 Let W be a < C px ' 1 -valued Schur function and let (p(z) = £nZ- Then the 
operator of composition by <p is a contraction from Ji?(W) into itself if and only if 
there exists a <C jXq -valued Schur function X(z) such that 



W(z)=W{e N z)X(z). 



(36) 
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Proof: By Theorem |l!6l the map Tq, is a contraction if and only if the kernel 

W( £jv z)W( £w w)*-W(z)W(w)* 

K w {z,w) ~K w (e N z,e N w) = 

1 — zw* 

is positive definite in the open unit disk. By Leech's factorization theorem, see loTl 
p. 107], the above kernel is positive definite if and only if there is a Schur function 
X{z) such that d36*b is in force. □ 

As an example, take any Schur function s and build 

W(z) = -±= (s(z) s(e N z) ■ ■ ■ sieg-h)) ■ (37) 

Then 

W(z) = W(e n z)Pn, 

where Pn is defined by (fJJ. 

7 The family tf N 

An effective approach to generating wavelet bases is the use of Multiresolution 
Analysis (MRA), see for example fPH [181 l26*ll . Traditionally one looks for a fi- 
nite family of functions in L2(R,c/x), or 1,2 (M. d ,dx) for some dimension d. If d = 1, 
one chooses a scale number, say N. If d > 1, instead one scales with ad x d matrix 
A over the integers. We assume that A is expansive, i.e., with eigenvalues bigger 
than 1 in modulus. If A is given, let N be the absolute value of its determinant. To 
create MRA wavelets we need an initial finite family J? of N functions in L2(K), 
or L?(M rf )- One of the functions is called the scaling function (<j) in the discussion 
below). For the moment, we will set d = 1, but the outline below easily generalizes 
to d > 1. An MRA wavelet basis is a basis for L,2(R), or L,2(M. d ) which is gener- 
ated from the initial family & and two operations : one operation is scaling by the 
number (or the matrix A if d > 1), and the other is action by integer translates 
of functions. The special property for the finite family of functions & is that if the 
Af-scaling is applied each function l// in F the result is in the closed span of the in- 
teger translates of the scaling function (p. The corresponding coefficients are called 
masking coefficients. The reason for this is that the scaled functions represent re- 
finements, and they are computed from masking points in a refinement. The role of 
the functions mo, mi,... ,mjv-i are the frequency response functions corresponding 
to the system of masking coefficients. From these functions we then build a matrix 
valued function W(z) as in d38l l. The question we address here is the characteriza- 
tion of the matrix valued function which arise this way. Now the wavelet filters we 
consider here go beyond those studied earlier in that we allow for wider families 
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of Multiresolution Analyses (MRAs). This includes more general wavelet families, 
allowing for example for wavelet frame bases, see e.g., |T3l HTl [42], multi-scale 
systems in dynamics, and in analysis of fractals; see (3TJ . 



7.1 The family characterization 

The filters we consider are matrix-valued (or operator valued) functions of a com- 
plex variable. In general if a positive integer A is given, and if a matrix function 
W(z) is designed to take values in C NxN , then of course, there are A 2 scalar-valued 
function occurring as matrix entries. However, in the case of filters arising in appli- 
cations involving A distinct frequency-bands, for example in wavelet constructions 
with scale number A, then we can take advantage of an additional symmetry for the 
given matrix function W(z), see for example (Q~|i in the Introduction. Here we point 
out that this A-symmetry condition (or A-periodicity) means that W(z) is then in 
fact determined by only A scalar valued functions, see ( 1381 below. These functions 
play three distinct roles as follows: They are (i) the scalar valued filter functions, Sj, 
for i = 0, 1, . . . ,N — 1, in generalized quadrature-mirror filter systems (the quadra- 
ture case corresponds to A = 2); they are (ii) scaling filters for scale-number A with 
each of the A scalar functions Sj generating an element in a wavelet system of func- 
tions on the real line and corresponding to scale-number A; and (iii) the system of 
scalar functions (si);=o,...jv-i generates an operator family (S;);=o,...jv-i constituting 
a representation of the Cuntz relations; thus generalizing Theorem 5.3 above. The 
results presented in this section are related to J9)- 

Recall that e N = «t. We shall say that a C A ' xA '-valued (A > 2) function W mero- 
morphic in the open unit disk D belongs to if it is of the form 



where sq, . . . , i are complex- valued functions meromorphic in D. Note that such 
a function, when analytic at the origin, will never be invertible there. A special case 
of this analyticity restriction of course is when W(z) has polynomial entries. Un- 
der the filter-to-wavelet[§correspondence lfl8l . polynomial filters are the compactly 
supported wavelets. In the sequel, it will turn out that we shall concentrate on the 
opposite cases. Namely, not only W(z) will not be analytic at the origin, in fact we 
shall have W(z)7 = Onxn- 

Recall that we have denoted by Pn the permutation matrix, 




/ S Q (z) S (e N z) 
s\ (z) Si (e N z) 



(38) 



\%-l(z) ^Ar-l(e^z) 



3 This correspondence: polynomial filter to compactly supported wavelet even works if d > 1 . 
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In-\ 0(JV— l)xl 
(see ©). 

Lemma 7.1 A < C NxN -valued function meromorphic in the open unit disk is of the 
form d38t if and only if it satisfies 

W(e N z) = W(z)P N 



Proof: Let W be a C A ' xA '-valued function meromorphic in B, and satisfying ([T), and 
let s\, . . . sn denote its columns, i.e. 



W(z) = (ji(z) j 2 (z) •••*iv(z)) 



(39) 



Namely, from 



/ *o(4 L z) \ 



1 



7 = 1, 



TV. 



Multiplying W by fly from the right makes a cyclic shift of the columns to the left, 
namely 

W(z)P N = (s 2 (z) s 3 {z) ■■■ s N {z) si(z)). 
Equation ((TJ then leads to 

(si(eivz) s 2 {e n z) ■■■ s N -i(e N z) s N (e N z)) = 

= (s 2 (z) S 3 (z) ••• Stf(z) Si(z)) . 

Thus 

*2(z) =ii(£ivz), s 3 (z) =si(e^z),...,s N (z) =si{e^ l z), 
and so W is of the asserted form. The converse is clear. □ 



Note that in contrast to Lemma I7TT1 in equation (l37l i we did not assume that W is 
square. 

When one assumes that the function W in the previous lemma is a generalized 
Schur function, the symmetry condition (fTJ can be translated into the realization. We 
present the result for the closely outer connected coisometric realization, but similar 
results hold for the closely inner connected isometric realization and connected uni- 
tary realizations as well (see Section R~2l for these notions). In the statement, recall 
that the state space & 1 will in general be infinite dimensional and endowed with a 
Pontryagin space structure. 
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Theorem 7.2 Let W be a generalized Schur function, and let 

W{z)=D + zC{I~zA)- l B 

be a closely-inner coisometric realization ofW, with state space SP. Then, W sat- 
isfies (HJ if and only if there is a bounded invertible operator T from into itself 
such that 



e N A B\ (T 0\ _ (T 0\ (AB 
e N CD lo/J ~ [o I N [CD 



(40) 



Furthermore, the operator T satisfies: 



T N =I. (41) 

Proof: The first equation follows from the uniqueness of the closely connected 
coisometric realization. Iterating d40b and taking into account that eft = 1 we get 

A B\ (T N 0\ _ (T N 0\ (AB 
CDj{o In) ~ \ In) \C D 

By uniqueness of the similarity operator we have T N = I. □ 



Proposition 7.3 Let W\ and W2 be in ^n- Then the functions 

W 1 {z)W 2 {z)* and Wi{z)Wz{\/z)* 
are meromorphic functions ofz N - 

Proof: Let W(z) = W\ (z)W 2 (z)*. Since P N P* = I N , we have 

W(e N z)=W l (e N z)W 2 {miy 
= W 1 (z)P N P N W 2 (z)* 
= W 1 (z)W 2 (z)* 
= W(z), 

that is 

W(e NZ ) = W(z). (42) 

The function W\ and W 2 are meromorphic in the open unit disk, and so is the function 
W. We denote by A the set of poles of W and by An the set of points w in the open 
unit disk such that vtA € A. Let for z = re' e with r > and 6 6 (— n, n], 

The function R is analytic in D \ {An U (— 1,0]}. Thanks to d42l it is continuous 
across the negative axis at those points in (—1,0) which are not in D\Ajy. It follows 
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that R is analytic in D \ An U {0}. Furthermore, W(z) = R(z N )- Any singular point 
of R is a pole (otherwise its roots of order N would be essential singularities of W), 
and so R is meromorphic in ID). □ 



In the rational case, the previous result has an easier and more precise proof. Indeed 
consider the Laurent expansion at the origin of W: 

W(z) = £ W k z k . 
-mo 

It converges in a punctured disk < |z| < r for some r > 0. Equation (142b implies 
that 

-m -m 

By uniqueness of the Laurent expansion we get that 



W k = 0, for k <£ NZ. 

Thus, if m > 0, we may assume without loss of generality that mo = Nhq for some 
no £ N. The function 

-l 

W-(z)= £ W k z k 

k=-m Q 

is rational, and so is the function 

W+(Z) = 

k=0 

We see that 

HL(z)= W "NZ" N 

—ma<nN<—N 

and so W_ (z) = /?_ (z^), where the function 

R-(z)= £ W„ w z" 

—mo<nN<—N 

is rational and analytic at infinity. The function W + is analytic at the origin, and thus 
can be written in realized form as: 

W+(z) =D + zC(I p -zA)- l B. 

Comparing with 



W+(z)= £\W 

n=0 



we have that 
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CA P B 



if p + 1 £ NN, 

W„ N if p + 1 = n2V, n E N. 



It follows that W + (z) = R+(^), where R + is the rational function defined by 



R+(z)=D+Y J z"CA 



rir^ \ nN— 1 



B 



n=l 



{n-l)N,N-l, 



The function 



is rational. 



= D + zC(I p -zA 
R(z)=R-(z)+R + (z) 



N '- l A N - l B. 



The proof of the preceding proposition can be mimicked to obtain the following 
result: 

Proposition 7.4 Let W\ and W2 be in *^v, with non identically vanishing determi- 
nant. Then there exists a me romorphic function R such that 



W i (z)W 2 (z)- l ^R(z N ). 
To this end, recall that the unitary matrix fjy, 



(43) 



Fn := ~7Tr 



( £- (0 - 0) 

c -(i-o) 

e -(2'0) 



-(0-1) 
-(hi) 
-(2-1) 



-(0-2) 
-(1-2) 
-(2-2) 



p -((AC-l)-0) P -((W-1)-1) -((iV-l)-2) 



e -(0-(JV-l)) ^ 
g -(l-(JV-l) 
e -(2-(JV-l) 

-((N-l)-(N-l) J 



generates the discrete Fourier transform. Namely, the discrete Fourier transform of 
x € is given by X = F^x and the inverse discrete Fourier transform is given by 
x = F£jX. Let furthermore 

W N (z) :=diag{l, z~\ ... , z l - N }F N . (44) 

With this special choice of W2 the previous proposition becomes: 
Proposition 7.5 W £ and detW ^ if and only if it can be written as 

W{z)=R(z N )W N {z) 1 

where R and Wn are as in ( |43t and ( 1441 ) respectively. 
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Let 



so that 



^( Z )=diag(^,/- 1 ^- 1 ,/- 2 4- 2 ,. 



■ ,ZEn), 



D N (l)=dmg(l,e%-\e k N 2 ,...,e N ). 



Functions which satisfy the related symmetry 

W(e NZ )^D N (iy 1 W(z)P N 



(45) 



appear in the theory of periodic systems. A function W satisfies ((45) if and only if 
it is of the form 



si (z) 



W(z) = 



So{e N z) 
si (e N z) 



%-l(z) -JrYSN-iiENZ) 



— js N -i(e N z) 



(46) 



See J5j Theorem 4. 1 , p. 38 1 ] . We note that the corresponding general bitangential in- 
terpolation problem (see flT4l for references) was solved in (5) for functions analytic 
and contractive in the open unit disk (that is, for Schur functions). Let us denote by 
SPer^ the family of functions meromorphic in the open unit disk and which satisfy 

gU). 

Proposition 7.6 The map W H> D^W is one-to-one from 3^er^ onto ^\r. If W is 
analytic and contractive in the open unit disk so is D^W . 

Proof: We first note that 

D N (e N z)=D N (z)D N (l). (47) 

Let now W € S^er^. In view of d47| i and ( |43T > we have 

D N (e N z) W (e N z) = D N (z)D N ( 1 )D N ( 1 ) - 1 W (z)P N 
= D N (z)W(z)P N , 

and so D N W £ c <o N . □ 

Epilogue: A reason for the recent success of wavelet algorithms is a coming to- 
gether of tools from engineering and harmonic analysis. While wavelets now enter 
into a multitude of applications from analysis and probability, it was the incorpo- 
ration of ideas from signal processing that offered new and easy-to-use algorithms, 
and hence wavelets are now used in both discrete problems, as well as in harmonic 
analysis decompositions. Following this philosophy we here employed tools from 
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system theory to wavelet problems and tried to show how ideas from wavelet de- 
compositions throw light on factorizations used by engineers. 

Since workers in wavelet theory often are not familiar with filterers in general, and 
FIR filters (short for Finite Impulse Response) in particular, widely used in the en- 
gineering literature, we have taken the opportunity to include a section for mathe- 
maticians about filters. Conversely (in the other direction), engineers are often not 
familiar with wavelet analysis, and we have included a brief exposition of wavelet 
facts addressed to engineers . We showed that there are explicit actions of infinite- 
dimensional Lie groups which accounts for all the wavelet filters; as well as for other 
classes of filters used in systems theory. Moreover, we described these groups, and 
explained how they arise in systems. The corresponding algorithms, including the 
discrete wavelet algorithms are used in a variety of multi-scale problems, as used for 
example in data mining. These are the discrete algorithms, and we described their 
counterparts in harmonic analysis in standard L2 Lebesgue spaces, as well as in re- 
producing kernels Hilbert spaces. We also outlined the role of Pontryagin spaces in 
the study of stability questions. 

In the engineering literature the study of filters is mostly confined to FIR filters. Re- 
call that FIR filters correspond to having the spectrum at the origin. In our previous 
work 1 9] we have explained that the set of FIR wavelet filters is small in a sense we 
made precise. This suggests two possible conclusions, 

(i) It is unrealistic to offer optimization schemes, over all FIR wavelet filters as part 
of the design procedure. 

(ii) It calls upon using, at least in some circumstances, also stable IIR (short for infi- 
nite impulse response) wavelet filters, i.e. the spectrum is confined to the open unit 
disk. 

The above extension to °^ 1n allows us to consider filters whose spectrum is in C \ T. 
The generalization to permits the spectrum to be everywhere in the complex 
plane. 

Roughly, we hope that this message will be useful to practitioners in their use of 
these rigorous mathematics tools. We offer algorithms hopefully improving on those 
used before. 
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